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1  Introduction 


It  now  is  recognized  widely  that  one  of  the  most  convenient  means  of  predicting 
the  unsteady  physics  of  turbulent  reacting  flows  is  via  large  eddy  simulation  (LES) 
[25,  37,  67].  A  serious  issue  in  such  simulations  is  accurate  modeling  of  the  subgrid 
scale  (SGS)  quantities.  While  this  issue  is  important  in  any  LES,  it  is  particularly 
difficult  when  dealing  with  reacting  flows  [37,  67,  5]. 

The  predictive  capability  of  probability  density  function  (PDF)  methods  in  Reynolds 
averaged  simulation  (RAS)  is  established  very  well  [67, 19,  32],  This  capability  is  due 
to  the  inherent  property  of  the  PDF  to  include  the  complete  statistical  information 
about  its  variables.  It  appears  that  the  use  of  PDF  methods  for  SGS  modeling  was 
suggested  first  by  Givi  [24];  however,  it  was  the  formal  definition  of  the  “filtered 
density  function”  (FDF)  by  Pope  [65]  that  facilitated  the  implementation  of  LES 
via  PDF.  One  of  the  first  implementations  of  LES/FDF  was  reported  by  Madnia 
and  Givi  [47]  in  which  the  shape  of  the  FDF  was  specified  a  priori.  This  “presumed 
FDF”  procedure  was  pursued  in  several  subsequent  studies,  in  most  of  which  it  was 
assumed  that  the  thermo-chemical  variables  depend  on  only  the  mixture  fraction, 
e.g.  infinitely  fast  reaction,  equilibrium  chemistry,  etc.  Therefore,  the  FDF  is 
univariate  [47,  11,  71,  39,  15,  16,  18,  40].  For  LES  of  non-equilibrium  reactive  flows, 
it  is  necessary  to  assume  the  “joint”  FDF  of  multi-scalars  [20,  21].  All  of  the  assumed 
SGS  scalar  FDFs  in  these  contributions  are  based  on  the  first  and  the  second  order 
moments.  The  PDFs  generated  in  this  way  facilitate  affordable  LES.  However,  it 
now  is  well  understood  that  the  true  PDF  depends  strongly  on  the  actual  physics 
of  mixing  in  a  given  flow  condition  [35]. 


1.1  Our  Previous  Work  in  LES/FDF 

A  reliable  means  of  determining  the  FDF  is  via  solution  of  its  transport  equation.  A 
systematic  study  in  this  endeavor  was  initiated  by  Colucci  et  al.  [10],  who  developed 
and  solved  a  modelled  transport  equation  for  the  marginal  scalar  FDF  (SFDF)  in 
constant  density  reacting  flows.  This  work  demonstrated,  for  the  first  time,  the 
feasibility  of  LES/FDF  and  the  importance  of  SGS  scalar  fluctuation  for  accurate 
modeling  of  the  filtered  reaction  rate.  The  extension  of  this  work  was  conducted 
by  Jaberi  et  al.  [34,  36],  who  developed  the  marginal  scalar  filtered  mass  density 
function  (SFMDF).  The  SFMDF  is  essentially  the  mass-weighted  form  of  the  SFDF. 
Due  to  the  marginal  nature  of  the  FDFs  in  these  studies,  all  of  the  hydrodynamic 
effects,  including  all  of  the  velocity-scalar  correlations,  are  modelled  via  conventional 
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(non-FDF)  methods.  The  extension  for  modeling  of  the  velocity  field  was  done  by 
Gicquel  et  al.  [23],  who  developed  and  solved  a  transport  equation  for  the  marginal 
velocity  FDF  (VFDF)  in  constant  density  flows.  This  work  demonstrated  some  of 
the  advantages  of  the  FDF  in  comparison  to  conventional  methods  in  accounting 
for  the  effects  of  SGS  velocity  correlations. 

1.2  Related  Work  by  Others 

Since  its  original  conception  [24,  65],  the  FDF  has  become  very  popular  for 
combustion  applications.  In  addition  to  our  previous  work,  as  cited  above,  the 
methodology  has  experienced  widespread  usage  by  many  others.  Examples  are 
contributions  in  its  basic  implementations  [90,  26,  70],  fine-tuning  of  its  sub-closures 
[72,  8,  31],  and  its  validation  via  laboratory  experiments  [82,  86,  83,  87].  The 
methodology  is  finding  its  way  into  commercial  codes  [26,  9]  and  has  been  the 
subject  of  detailed  discussions  in  several  recent  textbooks  [67,  19,  32], 


2  Accomplishments 

The  goal  of  this  research  was  to  improve  the  capabilities  of  the  FDF  method  and 
to  implement  it  for  LES  of  chemically  reacting  turbulent  flows.  We  feel  that  we 
have  been  very  successful  in  achieving  the  specific  objectives  of  this  work.  These 
objectives  were: 


1.  development  and  implementation  of  the  joint  velocity-scalar  filtered  density 
function  (VSFDF)  for  SGS  modeling  of  turbulent  combustion, 

2.  fine-tuning  and  implementation  of  the  scalar  FDF  (SFDF)  for  LES  of  complex 
flames  and  comparison  with  experimental  data. 


The  efforts  pertaining  to  these  objectives  are  described  partly  in  Refs.  [79]  and 
[78]  and  discussed  further  below.  In  addition,  the  PI  has  provided  an  invited  survey 
article  [25]  with  a  review  of  FDF  methods.  The  discussions  provide  a  comprehensive 
description  of  our  accomplishments  but  do  not  emphasize  the  materials  published 
in  articles  [79]  and  [78]. 
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2.1  Fundamental  Formulation 


In  a  compressible  flow  undergoing  chemical  reaction  involving  Ns  species,  the 
primary  transport  variables  are  the  density  p(x,  i),  the  velocity  vector  Uj(x,£) 
(i  =  1,2,3),  the  pressure  p(x,f),  the  total  specific  enthalpy  h(x,t)  and  the  species’ 
mass  fractions  ya(x,  f)  (a  =  1,  2, . . . ,  Ns).  The  equations  that  describe  the  transport 
of  these  variables  in  space  (xj)  (i  =  1,2,3)  and  time  ( t )  are  the  continuity, 
momentum,  enthalpy  (energy),  and  species’  mass  fraction  equations,  along  with 
an  equation  of  state 


dp  dpuj 

dt  +  dxj 

=  0, 

(la) 

dpui  dpujUi 

dt  '  dxj 

dp  drji 

dxi  dxj  ’ 

(lb) 

dp(j)Q  dpUjtpa 

dt  dxj 

dJa 

—  —  o  +  pSa,  a  —  1,2,  —  Ns  +  1, 

dxj 

(lc) 

V 

N, s 

=  pR0Tj2Yo/Ma  =  pRT, 

(Id) 

OL  —  1 


where  R0  and  R  are  the  universal  and  mixture  gas  constants  and  Ma  denotes  the 
molecular  weight  of  species  a.  The  chemical  reaction  source  terms  Sa  =  Sa  (<f>(x,t)) 
are  functions  of  compositional  scalars  (<f>  =  [4>i,4>2>  ■  •  •  ,4>ns+i])-  Equation  (lc) 
represents  the  transport  of  species’  mass  fraction  and  enthalpy  in  a  common  form 
with 

Ns 

4*a  —  Iqtj  Q!  Ns>  4*0  =  h  ^  (2) 

a=  1 


and 

h«  =  hl+  f  cr„(T')dT'.  (3) 

2  T0 

Here  T  and  To  denote  the  temperature  field  and  the  reference  temperature, 
respectively.  In  this  equation,  /i°  and  cPa  denote  the  absolute  enthalpy  at  To  and  the 
specific  heat  at  constant  pressure  for  species  a.  For  a  Newtonian  fluid,  with  Fick’s 
law  of  diffusion,  the  viscous  stress  tensor  and  the  scalar  flux  are  represented 
by 


Ja 


(  dui  duj  2  duk  .  \ 
\dxj  ^  dxi  3  dxk  v )  ’ 
dcj)Q 
1dxj ’ 


(4a) 

(4b) 
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<  * 

where  p  is  the  fluid  dynamic  viscosity  and  7  =  pT  denote  the  thermal  and  mass 
molecular  diffusivity  coefficients  for  all  scalars.  In  reactive  flows,  molecular  processes 
are  much  more  complicated  than  portrayed  by  Eq.  (4).  Since  the  molecular  diffusion 
typically  is  less  important  than  that  of  SGS  turbulence,  this  simple  model  is  adopted 
with  the  justifications  and  caveats  given  in  Refs.  [44,  64]  and  [3]. 

Large  eddy  simulation  involves  the  spatial  filtering  operation  [67,  63,  52,  22,  75] 

/+OO 

Q(x',t)G(x',x)dx',  (5) 

-OO 

where  G(x',x)  denotes  a  filter  function,  and  (Q(x,t))e  is  the  filtered  value  of  the 
transport  variable  Q(x,t).  In  variable-density  flows  it  is  convenient  to  use  the 
Favre-filtered  quantity  (Q(x,t))L  =  {pQ)t/  (p)t  We  consider  a  filter  function  that 
is  spatially  and  temporally  invariant  and  localized.  Thus,  G(x',x)  =  G(x'  —  x)  with 
the  properties  (2(x)  >0,  S-00  G(x)dx  =  1.  Applying  the  filtering  operation  to  Eqs. 
(1)  and  using  the  conventional  LES  approximation  for  the  diffusion  terms,  we  obtain 


d  ( p)(  (ui)  1  @_(p)e  (uj) L  (Ui)  L 


dt 


dx-j 


d  (p)e  (^0)4  _j_  d  (p)e  (uj)  1  ((fra)  1 


dt 


dxj 


d  (Ph 

dxi 

2d_ 

3  dxi 
_d_ 

dx 
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_d_ 

+  9xj  L 
d(uj) 

^  dx, 

d  {(Pa)  1 


P 


l  _|_  d  (uj)  1 


dxj 


dJVi) 

dxj  dxx 

_  d  {p)tTL[Uj,Uj) 
dxj 

d(p)irt(  Uj,4>0 

dx* 


(6a) 

+  (pSa)e  > 

(6b) 


where 


n{a,b)  =  (ab) L  -  (a) L  (b) L  (7) 

The  “velocity-scalar  filtered  mass  density  function”  (VSFMDF),  denoted  by 
Fi  ( v ,  ip;  x,  t),  is  defined  formally  as  [65] 

/  +  OO 

p(x',  t) c  {v,  </>;  u(x',t ),  cp(x',  t))  G(x'  -  x)dx',  (8) 

-OO 


where 


3  a 

(  (»>  w(x,  t),  0(x,  t))  =  JJ  S  (Vi  -  ut(x,  f))x[|(S  (ipa  -  <f>a(x,  <)) .  (9) 

2—1  a=l 

In  this  equation,  5  denotes  the  Dirac  delta  function,  and  v,  ip  are  the  velocity  vector 
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and  the  scalar  array  in  the  sample  space.  The  term  £  is  the  “fine-grained”  density 
[60,  64].  Equation  (8)  defines  the  velocity-scalar  filtered  mass  density  function 
(VSFMDF)  as  the  spatially  filtered  value  of  the  fine-grained  density.  With  the 
condition  of  a  positive  filter  kernel  [84],  FL  has  all  of  the  properties  of  a  mass  density 
function  [64].  For  further  developments  it  is  useful  to  define  the  “conditional  filtered 
value”  of  the  variable  Q(x,t)  as 


(Q(x,t)  u(x,t)  =  v,<f>(x,t)  -  ^  =  (Q  = 

Q  C*',  t)  p(x',  t) C  (v,  -0;  u(x*,  <),  <ft(x',  t))  G  (x*  -  x)  fix' 
FL(v,ip;x,t) 


(10) 


To  develop  the  VSFMDF  transport  equation,  we  consider  the  time  derivative  of  the 
fine-grained  density  function,  (Eq.  (9)).  Filtering  of  this  equation  yields  [77] 


dFL  dvtFL 
dt  dxi 


+ 


+ 


~wJSM)Fl] 

d_  //__!  dp 
dvi  \\p(4>)  dxi 

j}_(  (  1  ®l2i 

dvi  \\p(<j>)  dxj 

_d_  (/  1  dJl 
dlpa  \\p(<t>)  dxi 


V,l/> 


V,1p 


v,lp 


(11) 


which  is  an  exact  transport  equation  and  indicates  that  the  effects  of  convection 
(second  term  on  LHS)  and  chemical  reaction  (the  first  term  on  RHS)  appear  in 
closed  form.  The  unclosed  terms  denote  convective  effects  in  the  velocity-scalar 
sample  space.  Integrating  this  equation  over  the  v  field  yields  the  transport  equation 
for  the  scalar  filtered  mass  density  function  (SFMDF)  Fi(ip\x,t), 


dFjjjy,  x,  t) 
dt 


d[(rij(x,  t)|V>)<Fi,(V>;x,t)] 
dxi 


+ 


d 

d'lpa 
8  ' 
d-lpa 


[SQ(xl>)FL(-iP;x, 

/  1  dJl  \ 

\p{4>)  dxi  W / 


t)\ 

)  iv(Vs 


x,t) 

(12) 
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which  is  an  exact  transport  equation  for  the  SFMDF.  The  first  term  on  the  right 
hand-side  of  this  equation  is  due  to  chemical  reaction  and  is  in  closed  form.  The 
unclosed  nature  of  SGS  convection  and  mixing  is  indicated  by  the  conditional  filtered 
values. 

For  closure  of  the  VSFMDF  transport  equation,  following  previous  work  [28,  17, 
10,  23,  79],  we  utilize  the  generalized  Langevin  model  (GLM)  and  the  linear  mean 
square  estimation  (LMSE)  model  [60].  We  assume  a  constant  value  for  p,  =  7;  i.e. 
unity  Schmidt  (Sc)  and  Lewis  (Le)  numbers.  With  these  assumptions,  the  model 
is: 


dXf 


dW 


Uf  dt  + 


2  jJL 


( P ) 


dWi. 


1  d(p)L  |  2p  d2(uj)L  | 


P 


d 2  (uj)1 


(p)t  dxi  (p)t 


dxjdxj 


(p)e  dxjdxi 


(13a) 


2  P  d2(uj)L 

3  (p)t  dxidxj 


dt 


d<t> 


+ 

a 


+  gv  (uf  -  (Uj)L)  dt+  ^edw;  + 

=  -C^lo  (<p+ -  (<j)a)L)  dt  + Sa((f>+)dt, 


(13b) 

(13c) 


where  X+,  U+ ,  </>+  are  probabilistic  representations  of  position,  velocity  vector,  and 
scalar  variables,  respectively.  Also: 


G 


o 


U) 

e 

k 


1 

-OJ  I  r 


t 

k 

C' 


£3/2 

A,. 


:Cr 


6ij 


-jTL  ( Uk,Uk ) 


(14) 


Here  u>  is  the  SGS  mixing  frequency,  e  is  the  SGS  dissipation  rate,  k  is  the  SGS 
kinetic  energy,  and  A^  is  the  LES  filter  size.  The  parameters  Co,  C $  and  Ce  are 
model  constants  and  need  to  be  specified.  Thus,  the  Fokker-Planck  equation  [74]  is: 


1  a<p),8Ft 

(p)t  dxi  dvi  (p)[  dxjdxj  dvi  (p)t  dxidxj  dvi 

2  p  d^u^dFL  d  [(%  —  Fl\  &  (FJ_  (p)t) 

3  (p)[  dxjdxi  dv{  lJ  dvi  "  dx3dx3 

d_  (2 fid(ui)LdFL\  p  dJy^dW^  (¥Fl  lc  t  a2FL 

dxj  \(p)i  dxj  dvi )  (p),  dxk  dxk  dvidvj  2  (>t  dv3dv3 

c  J[{j>a-{4>a)L)FL\  _  djS^PW 

*  d'lpa  dlpa 

The  transport  equations  for  the  filtered  variables  are  obtained  by  integration  of  Eq. 
(15).  For  the  marginal  SFMDF,  the  model  for  the  scalar  transport  is  the  same.  But 
the  velocity  field  must  be  obtained  by  other  (non-FDF)  means.  In  this  case,  the 
physical  transport  is  modelled  via: 


OFl  dVrFL 
dt  dxi 

+ 

+ 

+ 


dXi(t)  = 


(Ui)l  + 


1  3(7  +  7*) 


(P)t 


dxi 


dt  +  ^2(7  +  'Yt)/{p)tdWi(t)  (16) 


which  yields  the  modelled  transport  equation  for  the  SFMDF: 


dFL  ,  d[(ut)LFL 


dt 


dxi 


d_ 

dxi 


(7  +  7t) 


djFJjph) 

dxi 


dip, 


dip0 


(17) 


Here,  7 1  is  the  SGS  diffusivity  and  must  be  modelled.  The  modified  kinetic  energy 
viscosity  (MKEV)  model  of  Jaberi  et  al.  [34]  has  proven  effective  for  this  purpose. 
The  model  in  this  way  is  equivalent  of  the  gradient-diffusion  model  for  the  SGS  flux 
of  the  scalars: 

d{(Pa)L 


(p)erL{uj,(pa)  =  -7 r 


dx3 


(18) 


The  mixing  frequency  here  is  denoted  by  Ctm  to  distinguish  it  from  that  in  the 
VSFMDF  because  it  has  to  be  modelled  differently.  The  model  proposed  by  Colucci 
et  al.  [10] 


Q>(7  +  7j) 

(p)AI 


(19) 


has  been  used  in  almost  all  of  the  previous  applications  of  SFMDF. 


From  the  computational  standpoint,  solution  of  the  SDEs  (e.g.  Eq.  (13)  and/or  Eq. 
(16)  is  significantly  easier  than  the  modelled  FDF  transport  equations  (15),  (17). 
The  most  effective  way  of  solving  Eq.  (13)  and/or  Eq.  (16))  is  via  Monte  Carlo 
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(MC)  methods.  These  methods  have  been  used  for  simulation  of  a  wide  variety  of 
stochastic  problems  [27]  and  have  benefited  significantly  from  modern  developments 
in  SDE  solver  technology  [42].  MC  methods  have  been  the  primary  means  of  solving 
the  PDF  in  RAS  [64,  29,  30,  66,  89]  and,  thus  far,  the  primary  method  of  choice 
for  solving  the  FDF  in  LES.  Typically,  the  method  is  implemented  by  representing 
the  FDF  by  an  ensemble  of,  say  Np,  particles.  These  particles  carry  information 
pertaining  to  their  positions,  X^(i),  velocities,  U^(t),  and  scalar  values,  t ), 

n  =  1, . . . ,  Np.  This  information  is  updated  via  temporal  integration  of  the  modelled 
SDEs.  While  it  is  potentially,  or  eventually,  possible  to  simulate  FDF  exclusively 
via  MC,  the  most  practical  procedure  is  via  “hybrid”  methods  in  conjunction  with 
“deterministic”  schemes.  There  has  been  significant  progress  in  the  development 
of  high-order  finite  difference  (FD),  finite  volume,  and  spectral  methods.  A  hybrid 
method  would  make  use  of  this  high-order  accuracy  which  is  not  yet  achievable 
in  methods  exclusively  via  MC.  Moreover,  the  influence  of  the  MC  dispersion  and 
statistical  errors  is  less  significant  in  hybrid  methods.  Finally,  in  many  cases  it  is 
easier  to  specify  the  model  parameters  in  a  deterministic  fashion  than  via  MC.  These 
issues  are  investigated  in  detail  in  the  context  of  RAS/PDF  [57,  59]  and  constitute 
a  major  element  of  the  computational  procedure  in  LES/FDF  [34,  23,  79]. 

The  base  flow  simulations  are  based  on  a  compact  parameter  finite-difference  (FD) 
scheme  [7,  41]  on  a  number  of  fixed  grid  points  with  spatial  spacing  A.  The 
MC  particles  are  distributed  randomly  and  are  free  to  move  anywhere  within  the 
domain.  The  statistical  information,  e.g.  filtered  values,  at  any  point  is  obtained  by 
considering  an  ensemble  of  NE  particles  residing  within  an  ensemble  domain  of  side 
length  A E  centered  around  the  points.  For  reliable  statistics  with  minimal  numerical 
dispersion,  it  is  desired  to  minimize  the  size  of  ensemble  domain  and  maximize  the 
number  of  the  MC  particles  [64].  In  this  way,  the  ensemble  statistics  would  tend 
to  the  desired  filtered  values.  Transfer  of  information  from  the  grid  points  to  the 
MC  particles  is  accomplished  via  interpolation.  The  transfer  of  information  from 
the  particles  to  the  grid  points  is  accomplished  via  ensemble  averaging  as  described 
above. 

The  FD  solver  determines  the  pressure  field  which  is  used  further  in  the  MC 
procedure.  The  transport  equations  solved  by  the  FD  include  unclosed  second 
order  moments  which  are  obtained  from  the  MC  solver.  The  FD  solver  also 
determines  the  filtered  velocity  and  scalar  fields.  That  is,  there  is  a  “redundancy” 
in  the  determination  of  the  first  filtered  moments  as  both  the  FD  and  the  MC 
procedures  provides  the  solution  of  this  field.  This  redundancy  actually  is  very 
useful  in  monitoring  the  accuracy  of  the  simulated  results  as  shown  in  previous 
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work  [34,  57,  59,  23,  79]. 


To  reduce  the  computational  cost,  a  procedure  involving  the  use  of  non-uniform 
weights  [34]  also  is  considered.  This  procedure  allows  a  smaller  number  of  particles 
in  regions  where  a  low  degree  of  variability  is  expected.  Conversely,  in  regions  of  high 
variability,  a  large  number  of  particles  is  allowed.  It  has  been  shown  [64,  34]  that  the 
sum  of  weights  within  the  ensemble  domain  is  related  to  filtered  fluid  density,  and 
the  Favre-filtered  values  are  constructed  from  the  weighted  average  values.  With 
uniform  weights  [64],  the  particle  number  density  decreases  significantly  in  regions 
of  low  density,  such  as  reaction  zone.  The  implementation  of  variable  weight  allows 
the  increase  in  particle  density  without  increasing  the  particle  number  density  in 
these  regions. 


2.2  Latest  Results  via  VSFMDF 

Simulations  are  conducted  of  a  three-dimensional  temporally  developing  mixing  layer 
involving  transport  of  a  passive  scalar.  The  VSFMDF  predictions  are  compared 
with  data  obtained  by  direct  numerical  simulation  (DNS)  of  the  same  flow.  For 
a  comparative  assessment  against  a  conventional  method,  LES  of  this  layer  also  is 
performed  using  the  Smagorinsky  [80]  SGS  closure.  In  the  representation  below, 
x,  y ,  and  z  denote  the  streamwise,  the  cross-stream,  and  the  spanwise  directions, 
respectively.  The  velocity  components  along  these  directions  are  denoted  by  u ,  v, 
and  w  in  the  x,  y,  and  z  directions,  respectively.  The  temporal  mixing  layer 
consists  of  two  parallel  streams  travelling  in  opposite  directions  with  the  same  speed 
[73,  76,  56].  The  normalized  filtered  streamwise  velocity,  scalar,  and  temperature 
fields  are  initialized  with  hyperbolic  tangent  profiles  with  (u)L  =  1,  ((f>)L  =  1, 
(T)l  =  1  for  the  top  stream  and  (u)L  =  —  1,  (<f>)L  =  0,  (T)L  =  2  for  the  bottom 
stream.  The  length  Lv  is  specified  such  that  Lv  =  2Np Xu,  where  Np  is  the  desired 
number  of  successive  vortex  pairings  and  \u  is  the  wavelength  of  the  most  unstable 
mode  corresponding  to  the  mean  streamwise  velocity  profile  imposed  at  the  initial 
time.  The  flow  variables  are  normalized  with  respect  to  the  half  initial  vorticity 
thickness,  Lr  =  }  where  { u)L  is  the  Reynolds-averaged 

value  of  the  filtered  streamwise  velocity  and  AU  is  the  velocity  difference  across  the 
layer). 

We  consider  a  cubic  box,  0  <  x  <  L,  ^  <  y  <  |,  0  <  z  <  L  where  L  =  Lv/Lr. 
The  3D  field  is  parameterized  in  a  procedure  somewhat  similar  to  that  in  Ref.  [85]. 
Simulations  are  conducted  on  equal  grid  spacings  Ax  =  Ay  =  Az  =  A  with  the 
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number  of  grid  points  1933  and  333  for  DNS  and  LES,  respectively.  The  number  of 
particles  per  grid  point  is  NPG  =  320  (Ne  =  40)  and  the  ensemble  domain  size  is 
set  equal  to  As  x  Ag  x  Ag  where  Ag  =  A/2.  For  comparison,  the  DNS  data  are 
filtered  from  the  original  high  resolution  to  the  coarse  points  using  a  top-hat  filter 
function  with  the  filter  size  (Ag)  equal  to  2A. 

The  formation  of  the  large  scale  structures  are  expedited  through  eigenfunctions 
based  initial  perturbations  [53,  45].  We  impose  two-dimensional  [76,  85,  54]  and 
three-dimensional  [76,  55]  perturbations  with  a  random  phase  shift  between  the 
3 D  modes.  This  forcing  results  in  the  formation  of  two  successive  vortex  pairings 
and  strong  three-dimensionality.  For  a  statistical  appraisal,  we  also  consider  the 
“resolved”  and  the  “total”  components  of  the  Reynolds-averaged  moments.  The 
former  are  denoted  by  R(a,b)  with  R(a,b)  =  ^(a)L  —  (a)LJ  ((b) L  —  ( b) L and  the 
latter  are  r(a,  b )  with  r(a,  b)  =  (a  —  a)  (b  —  6) .  In  DNS,  the  “total”  components  are 
available  directly,  while  in  LES  they  are  approximated  by  r(a,  b)  ~  R(a ,  b )  +  rL(a,  b) 
[85]. 

Figure  1  shows  the  instantaneous  iso-surface  of  the  (<j>)L  field  at  t  =  80  as  obtained 
by  DNS  and  by  VSFMDF  and  Smagorinsky  models.  By  this  time,  the  flow  is  going 
through  pairings  and  exhibits  strong  3D  effects  [53].  The  two  neighboring  rollers  are 
being  paired,  and  the  formation  of  secondary  structures  is  evident.  As  illustrated  in 
this  figure  and  consistent  with  the  previous  results  [23],  the  amount  of  SGS  diffusion 
with  the  Smagorinsky  model  is  significant.  Thus,  the  predicted  results  are  overly 
smooth,  while  there  is  more  resemblance  between  the  VSFMDF  and  DNS  results. 

Several  components  of  the  Reynolds-averaged  values  of  the  second  order  SGS 
moments  are  compared  with  DNS  data  in  Fig.  2.  In  general,  the  VSFMDF  results 
are  in  better  agreement  with  DNS  data  than  those  predicted  by  the  Smagorinsky 
model.  In  this  configuration,  there  are  no  strong  velocity  and  scalar  gradients  in  the 
strcamwise  and  spanwise  directions,  and,  thus,  a  gradient-diffusion  type  model  such 
as  Smagorinsky  is  not  capable  of  providing  correct  prediction  of  scalar  flux  values 
in  these  directions.  Consequently,  the  VSFMDF  is  expected  to  be  more  effective  for 
LES  of  reacting  flows,  provided  that  the  extent  of  SGS  mixing  is  influenced  heavily 
by  these  SGS  moments  [4,  62],  Several  components  of  the  resolved  second  order 
moments  are  presented  in  Fig.  3.  As  expected,  the  performance  of  the  Smagorinsky 
model  is  not  satisfactory,  as  it  does  not  predict  the  spread  and  peak  values  accurately. 
The  VSFMDF  provides  more  reasonable  predictions.  The  “total”  components  also 
yield  very  good  agreement  with  DNS  data,  as  shown  in  Fig.  4.  As  observed  in  Refs. 
[79,  23],  the  “total”  components  predicted  by  VSFMDF  are  almost  insensitive  to 
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the  model  parameters. 


2.3  Latest  Results  via  SFMDF 

LES  via  SFMDF  was  conducted  of  a  turbulent  bluff-body  stabilized  hydrogen- 
methane  jet  flame.  This  flame  has  been  studied  experimentally  by  the  Combustion 
Research  Facility  at  the  Sandia  National  Laboratories  and  the  Thermal  Research 
Group  at  the  University  of  Sydney  [49,  48,  50,  14,  13,  12,  2,  1].  Bluff-body  stabilized 
flames  have  been  studied  by  several  investigators  [51,  38,  58,  46,  68,  69,  70,  33,  88, 
43,  12].  These  flames  produce  complex  flow  patterns  characteristic  of  practical 
combustors,  and  are  important  for  industrial  applications. 

In  the  experiments,  a  variety  of  operating  conditions  are  considered.  The  flow 
schematic  is  shown  in  Fig.  5.  The  central  round  fuel  jet  is  surrounded  by  a  round 
bluff-body  and  air  coflow.  The  jet  diameter  Dj= 3.6  mm,  and  the  bluff-body 
diameter  Dg= 50  mm.  The  recirculation  zone  immediately  follows  the  bluff-body 
surface.  The  neck  zone  provides  a  controlled  region  where  the  turbulent  mixing  rate 
is  significant  and  flame  blow-off  (extinction)  can  occur  [12].  Further  downstream, 
the  flame  exhibits  jet  like  behavior.  Both  non-reacting  and  reacting  flows  are 
considered.  The  experiments  pertaining  to  the  former  are  distinguished  by  their 
corresponding  bulk  jet  velocities,  Uj,  while  those  of  the  latter  are  distinguished  by 
their  blow-off  parameter.  The  recirculation  zones  differ  in  length  from  ~  1.0 Dg 
for  the  non-reacting  cases  to  ~  1.6 Dq  for  the  reacting  cases.  The  reacting  case  is 
characterized  by  a  Reynolds  number  of  15, 800  based  on  fuel  jet  diameter  and  bulk  jet 
velocity.  The  bulk  jet  and  the  coflow  velocities  are  118  m/s  and  40  m/s,  respectively. 
This  condition  corresponds  to  50%  blow-off.  The  100%  blow-off  velocity  of  the  flame 
is  /7bo=235  m/s. 

The  most  complete  set  of  measurements  has  been  compiled  for  the  hydrogen- 
methane  flames.  These  flames  are  considered  at  50%,  75%,  and  91%  blow-off 
(extinction).  The  jet  is  maintained  at  a  temperature  of  298  K  and  is  composed 
of  50%  hydrogen  (H2)  and  50%  methane  (CH4)  by  volume  [13].  The  dilution  of 
hydrogen  in  methane  reduces  the  formation  of  soot.  The  coflow  is  maintained  at 
room  temperature  and  is  composed  of  air.  The  bluff-body  is  a  ceramic  surface  that 
heated  up  to  an  average  of  1003  K  during  operation.  The  effects  of  radiative  heat 
transfer  on  the  temperature  and  other  major  species  are  small  [33]. 

LES  is  conducted  of  the  50%  blow-off  case.  The  chemistry  is  near-equilibrium  and 
is  parameterized  by  the  flamelet  library  constructed  by  simulation  of  a  counterflow 
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(opposed  jet)  laminar  flame  [61].  This  flame  is  modelled  with  detailed  chemical 
kinetics  [81,  6].  At  low  strain  rates,  x>  all  °f  the  thermo-chemical  variables  are 
related  to  the  “mixture  fraction”  Z(x,i);  that  is,  Q(x,t)  =  Q(Z(x,t)).  Therefore, 
(p(x,  t))f  (Q(x,  t))L  =  Q{'iI>z)F'l('iPz]  x,  t)  dipz-  The  flamelet  table  constructed 
in  this  way  is  used  in  conjunction  with  SFMDF  predictions  of  the  mixture  fraction. 

The  simulations  are  conducted  on  a  3D  Cartesian  mesh  with  uniformly  spaced  grid 
points.  The  computational  domain  spans  a  region  of  30D j  x  22.5 D j  x  22.5 D j  in  the 
streamwise  (x),  and  the  two  lateral  (y,  z)  directions,  respectively.  The  number  of 
grid  points  is  101  x  151  x  151  in  the  x,  y,  and  2  directions,  respectively.  The  filter  size 
is  set  at  =  2(AxAj/Az)P3,  where  the  Ax,  Ay,  and  Az  are  the  grid  spacings  in 
the  corresponding  directions.  The  size  of  the  ensemble  domain  in  the  MC  simulation 
is  equal  to  the  filter  size.  There  are  approximately  40  particles  in  each  ensemble 
domain.  Per  results  of  extensive  previous  studies  [34,  10,  23,  79],  this  number  is 
sufficient  to  yield  an  excellent  statistical  accuracy,  with  minimal  dispersive  errors. 
There  are  about  6.8  million  MC  particles  within  the  computational  domain  at  all 
times.  The  results  are  monitored  to  ensure  that  the  particles  fully  encompass  and 
extend  well  beyond  the  regions  of  non-zero  vorticity  and  reaction. 

The  statistics  are  obtained  by  long-time  averaging  of  the  Favre  filtered  fields  during 
several  residence  times.  Data  collection  is  initialized  after  the  flow  has  swept  the 
domain  during  the  initial  four  flow  through  times.  The  PDFs  are  constructed  based 
on  the  SFMDF  data  resolved  at  the  FD  grid  points.  The  notations  Q  and  RMS(Q) 
denote  the  time-averaged  mean  and  root  mean  square  values  of  the  variable  Q , 
respectively.  The  radial  direction  is  denoted  by  r  =  \/y^+  22.  The  flamelet  table 
with  a  constant  strain  rate  of  x  =  100  1/s  is  used.  This  choice  is  consistent  with 
previous  studies  of  bluff-body  flames  [58]  and  experimental  assessments  [13]. 

Tables  1  and  2  show  the  specifications  of  the  experiments  and  simulation  data, 
respectively.  The  experimental  flow  fields  of  the  reacting  cases  HMlE-S(l-2) 
considered  at  Sydney  are  meant  to  represent  the  same  flame  conditions  as  the 
HM1  experiments  considered  at  Sandia.  However,  the  wind  tunnel  in  Sydney 
did  not  provide  the  same  hydrodynamic  conditions.  The  streamwise  velocities  are 
Uj  —  118m/s  and  Uc  =  40m/s  in  the  HMl,  and  Uj  =  108 m/s  and  Uc  =  33m/ s 
in  the  HM1E-S  experiments.  The  new  jet  and  coflow  velocities  were  chosen  such 
that  the  HMl  and  HM1E-S  flames  are  equally  proportional  (within  50%  blow-off) 
to  their  corresponding  blow-off  velocities.  To  facilitate  comparisons  of  the  flow 
statistics  of  the  HM1E-S  data  with  those  of  HMl  and  SFMDF,  the  HM1E-S  values 
of  the  streamwise  velocity  are  scaled  by  a  ratio  of  the  bulk  jet  velocities;  that  is 
118/108. 
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Table  1:  Experimental  specifications  [1]. 


Label 

Fuel(RaticF) 

Uj* 

UCoc 

Available  data 

Year 

B4C1-S1 

Air 

61 

20 

U,  V , 

RMS{U ),  RMS(V) 

1995 

B4C1-S2 

Air 

61 

20 

U,  V,  RMS(U), 
RMS(V),  RMS(UV ) 

12/01/1999 

B4C1-S3 

Air 

61 

20 

U,  V,  RMS(U) 
RMS(V),  RMS(UV) 

12/15/1999 

B4C2-50 

CNG 

50 

20 

Z,  RMS{Z) 

B4C2-85 

CNG 

85 

20 

Z,  RMS(Z) 

B4C2-143 

CNG 

143 

20 

Z,  RMS(Z) 

HM1 

H2 :  CNG 
(1:1) 

118 

40 

u,  V, 

RMS(U),  RMS(V ) 

1995 

HM1 

>* — ^  to 
i — 1 

^  o 

w  Stj 

118 

40 

Z,  T,  YJ,  RMS(Z), 
RMS(T),  RMS(Ya), 
PDFs 

1995 

HM1E-S1 

H2  :  CNG 
(1:1) 

108 

35 

U,  V,  RMS(U), 
RMS{V),  RMS(UV) 

01/21/2000 

HM1E-S2 

H2  :  CNG 
(1:1) 

108 

35 

U,  V,  RMS{U ), 
RMS(V),  RMS(UV) 

02/11/2000 

t  Volumetric. 

{[m/s] 

§02,  N2,  Ht,  H20 ,  CO,  C02,  OH,  NO 


Table  2:  Simulation  specifications. 


Label 

Fuel(Ratio^) 

Uj* 

UCoc 

Re5  l 

i 

LES-61 

Air 

61 

20 

13899 

LES-50 

CHa 

50 

20 

10580 

LES-85 

CHa 

85 

20 

17987 

LES-143 

ch4 

143 

20 

30260 

SFMDF 

H2  :  CHa 
(1:1) 

118 

40 

15800 

*  Volumetric, 
tfm/s] 

§  Reynolds  number  based  on  Dj  and  Uj 
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The  non-reacting  experiments,  B4Cl-S(l-3),  consider  air  at  both  the  fuel  and  coflow 
streams.  The  bulk  streamwise  velocity  of  61  m/s  and  the  coflow  streamwise  velocity 
of  20  m/s  are  considered.  The  recirculation  zone  with  the  double  vortex  structures 
captured  by  the  experiments  and  LES  are  shown  via  the  time-averaged  streamwise 
velocity  contours  with  superimposed  streamlines  and  velocity  vectors  in  Fig.  6.  The 
model  captures  the  two  recirculation  zones  well.  The  length  of  the  recirculation 
region  behind  the  bluff-body  surface  is  overpredicted  slightly.  The  counter-rotating 
vortices  also  are  predicted  slightly  further  downstream  from  their  experimentally 
reported  counterparts. 

The  non-reacting  experiments,  B4C2,  consider  a  “cold”  methane  fuel  stream  with 
bulk  streamwise  velocities  of  50,  85,  and  143  m/s,  and  the  air  coflow  stream  with 
the  streamwise  velocity  of  20  m/s.  The  time  averaged  streamwise  velocity  with 
superimposed  streamlines  and  velocity  vectors  are  shown  in  Fig.  7.  As  observed 
previously,  a  set  of  two  counter-rotating  vortices  forms  behind  the  bluff-body  surface. 
The  recirculation  region  shortens  as  the  fuel  stream  velocity  increases.  In  addition, 
the  two  vortices  reposition  themselves,  with  the  inner  vortex  shifting  downstream 
with  respect  to  the  outer  one.  The  shift  is  accompanied  by  a  loss  in  inner  vortex 
circulation  strength  that  becomes  most  noticeable  in  LES-143.  In  the  latter  case 
large  streamwise  velocity  causes  the  inner  vortex  to  break  down  partially.  The 
instantaneous  and  the  time-averaged  3D  contours  of  the  mixture  fraction  are  shown 
in  Fig.  8.  Different  mixing  structures  are  observed  in  the  instantaneous  values  as  the 
streamwise  velocity  of  the  fuel  stream  increases.  The  time-averaged  plots  accentuate 
the  shortening  of  the  recirculation  zone  as  highlighted  by  the  “cone”  of  the  mixture 
fraction  iso-surface  of  0.06. 

The  capability  of  the  SFMDF  to  predict  the  reacting  flow  field  now  is  demonstrated. 
The  time  averaged  streamwise  velocity  contours  with  the  streamlines  and  velocity 
vectors  superimposed  are  shown  in  Fig.  9.  This  figure  shows  the  recirculation  region 
and  the  two  characteristic  counter-rotating  vortices.  The  recirculating  zones  are 
captured  well  by  the  simulations.  However,  there  are  a  few  discrepancies  that  are 
manifested  by  the  relative  position  of  the  inner  and  outer  vortices.  In  the  HMl 
experiment,  as  well  as  in  the  simulation,  the  inner  vortex  is  located  behind  the 
outer  vortex.  However,  in  the  HM1E-S  experiments  the  situation  is  reversed.  The 
recirculating  region  extends  to  ~  1.6 DB,  as  compared  to  ~  1 DB  in  the  isothermal 
cases.  The  central  jet  penetrates  more  easily  the  low  density  region.  The  mean 
profiles  of  the  streamwise  velocity  and  the  mixture  fraction  are  shown  in  Figs.  10 
and  11,  respectively.  The  comparisons  show  good  overall  qualitative  agreement  with 
data. 
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The  radial  profiles  of  the  mean  and  the  RMS  values  of  the  temperature  are  shown 
in  Figs.  12  and  13,  respectively.  The  near  field  mean  values  compare  well  with 
the  experimental  data.  The  downstream  locations,  however,  are  overpredicted  on 
the  oxidizer  stream  side  of  the  bluff-body.  The  overprediction  can  be  traced  to  the 
temperature  on  the  fuel  lean  side  of  the  flamelet  table.  The  mean  temperature  is  the 
highest  in  the  central  mixing  layer  region  where  the  hot  products  of  combustion  are 
convected  toward  the  bluff-body  by  the  counter-rotating  vortical  structures.  The 
experimental  temperature  RMS  values  are  overpredicted  by  the  total  RMS  values 
and  generally  are  predicted  well  by  the  resolved  RMS  values  with  the  exception 
of  the  oxidizer  stream  side  of  the  bluff-body  at  downstream  locations.  Similar 
to  the  mixture  fraction  RMS,  the  effects  of  the  three  mixing  layers,  generated 
by  the  counter-rotating  vortical  structures,  are  characterized  at  x/Db  —  0.8  and 
t/Rb  ~  0.2,  0.5,  0.9  by  the  peaking  RMS  values  of  the  temperature.  The  PDFs 
of  the  resolved  scalar  fields,  as  predicted  by  SFMDF,  are  compared  with  those 
measured  experimentally  in  Fig.  14  for  the  resolved  mixture  fraction  and  Fig.  15  for 
the  resolved  temperature.  The  overall  agreement  between  SFMDF  predictions  and 
the  experiment  is  very  good. 


2.4  Computational  Requirements 

For  the  computational  requirements  of  FDF,  the  relative,  normalized  run  times  for 
conventional  (non-FDF):  SFDF  :  VSFDF  :  DNS  are:  (1)  :  (3  -  5)  :  (15  -  20) 
:  (1500  —  2000).  Obviously,  the  most  serious  issue  associated  with  FDF  is  its 
computational  cost.  In  most  cases,  LES/FDF  is  more  expensive  than  conventional 
LES.  The  overhead  in  comparison  to  non-FDF  methods  is  expected,  considering 
all  of  the  SGS  statistical  information  that  LES/FDF  provides.  The  computational 
time  for  FDF  is  significantly  less  than  that  of  DNS.  It  is  not  claimed  that  LES/FDF 
is  an  alternative  to  DNS;  nor  is  it  claimed  that  FDF  is  capable  of  reproducing  all 
DNS  results,  but  this  comparison  could  be  done  for  cases  where  DNS  was  possible. 
The  close  proximity  of  values  obtained  by  LES/FDF  and  DNS  and  the  substantially 
lower  computational  costs  of  FDF  make  it  a  viable  tool  for  simulations  of  reacting 
flow  systems  for  which  DNS  is  not  possible.  With  the  development  of  more  efficient 
FDF  solvers,  it  is  predicted  that  LES/FDF  will  distinguish  itself  as  a  major  tool  for 
prediction  of  engineering  combustion  problems. 
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3  Interaction  with  AFRL 


We  are  very  grateful  to  have  the  opportunity  of  interacting  with  several  of  the 
researchers  at  the  Wright-Patterson  AFB  on  several  occasions.  The  last  time  such 
interactions  happened  was  in  the  summer  of  2005,  when  the  PI  visited  AFRL.  Per 
recommendation  of  Dr.  James  Gord  of  AFRL/PR,  the  PI  met  with  Drs.  Vish  Katta 
and  Terry  Meyer.  The  PI  also  met  Dr.  Datta  Gaitonde  of  AFRL/VA  and  gave  a 
tutorial/review  seminar  on  SGS  modeling  in  combustion  for  both  AFRL/VA  and 
AFRL/PR.  This  seminar  was  attended  by  more  that  50  of  AFRL  scientists.  It  is 
obvious  that  the  subject  of  LES  is  of  significant  interest  at  AFRL. 


4  Enhancement  of  Technology  and  Education 


With  completion  of  this  research,  we  have  been  able  to  contribute  to  maintain  U.S. 
leadership  in  a  technology  which  is  of  significant  importance  to  DOD  and  AFOSR.  In 
order  to  demonstrate  our  visibility  in  this  research,  here  we  shall  list  all  the  awards 
and  some  of  noticeable  achievements  of  the  personnel  involved  in  this  program. 


4.1  Graduate  Students 

Involvement  of  students  in  research  is  an  issue  which  is  taken  very  seriously  at  our 
University.  We  are  committed  to  recruiting  excellent  quality  students  and  involving 
them  in  high  caliber  research.  During  the  past  three  years,  the  following  students 
have  been  supported  under  this  Grant. 

1.  Tomasz  G.  Drozda,  Ph.D.  in  Mechanical  Engineering,  December  2005. 
Dissertation:  Implementation  of  LES/SFMDF  for  Prediction  of  Non- 
Premixed  Turbulent  Flame.  Current  Position:  Research  Scientist, 
Combustion  Research  Facility,  Sandia  National  Laboratories,,  Livermore, 
CA. 

2.  Reza  M.H.  Sheikhi,  Ph.D.  in  Mechanical  Engineering,  December  2005. 
Dissertation:  Joint  Velocity-Scalar  Filtered  Density  Function  for  Large 
Eddy  Simulation  of  Turbulent  Reacting  Flows.  Current  Position: 
Research  Assistant  Professor,  the  University  of  Pittsburgh,  Pittsburgh, 
PA. 


18 


4.2  Awards  and  Honors 


1.  Peyman  Givi  received  the  NASA  Public  Service  Medal:  For  contribution  in 
the  development  of  advanced  modeling  techniques  for  designing  high-speed 
and  hypersonic  propulsion  systems  for  aerospace  vehicles  (2005). 

2.  Peyman  Givi  was  named  Distinguished  Alumnus,  Phi  Kappa  Phi  Honor 
Society,  Youngstown  State  University,  Youngstown,  OH  (2004). 

3.  Peyman  Givi’s  Profile  was  featured  in: 

•  Pittsburgh  Tribune  Review ,  “Newsmaker,”  page  B.2,  September  8,  2005. 

•  Youngstown  State  University  Magazine,  “An  Outstanding  Graduate”, 
part  of  the  article  “60  Years  of  Engineering  Excellence  at  YSU,”  pp. 
22-23,  Summer/Fall  2004. 

•  Technology  Horizons:  “Aero  Propulsion  Combustors  -  A  New 

Modeling  Capability  Quantitatively  Predicts  Aero  Propulsion  Combustor 
Performance,”  5(1),  p.  46,  February  2004. 

•  Pittsburgh  Post- Gazette:  “Europe  Slow  in  Stemming  Brain  Drain  to 
America”  and  “Pittsburgh  Brain  Gain”,  77(81),  October  20,  2003.  Cover 
page  and  page  A7. 

•  Space  Directorate  Highlights,  publication  of  NASA  Glenn  Research 
Center:  Significant  Research  Finding:  Sub-Grid  Scale  Modeling  in 
Turbulent  Combustion,  June  21,  2003.  Cleveland,  OH. 


4.3  Publications 

In  all  of  the  following  publications,  the  support  from  AFOSR  is  gratefully 
acknowledged: 

Invited  Review  Articles  and  Book  Chapters: 


1.  P.  Givi,  “Filtered  Density  Function  for  Subgrid  Scale  Modeling  of  Turbulent 
Combustion,”  AIA A  Journal ,  44(1),  16-23  (2006). 

2.  P.  Givi,  “Subgrid  Scale  Modeling  in  Turbulent  Combustion:  A  Review,”  A1AA 
Paper  2003-5081 ,  2003. 
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3.  C.K.  Madnia,  F  A.  Jaberi  and  P.  Givi  “Large  Eddy  Simulation  of  Heat  and 
Mass  Transport  in  Turbulent  Flows,”  Chapter  5  in  Handbook  of  Numerical 
Heat  Transfer,  Second  Edition,  Editor:  W.J.  Minkowycz,  E.M.  Sparrow,  and 
J.Y.  Murthy,  John  Wiley  &  Sons,  Inc.,  New  York,  NY,  2005. 


Journal  Articles  in  Print: 


1.  M.R.H.  Sheikhi,  T.G.  Drozda,  P.  Givi,  F.A.  Jaberi  and  S.B.  Pope,  “Large  Eddy 
Simulation  of  a  Turbulent  Nonpremixed  Piloted  Methane  Jet  Flame  (Sandia 
Flame  D),”  Proceedings  of  the  Combustion  Institute,  30,  549-556,  2005. 

2.  M.R.H.  Sheikhi,  T.G.  Drozda,  P.  Givi  and  S.B.  Pope,  “Velocity-Scalar  Filtered 
Density  Function  for  Large  Eddy  Simulation  of  Turbulent  Flows,”  Physics  of 
Fluids ,  15(6),  2321-2337,  2003. 


Conference  Papers: 


1.  M.R.H.  Sheikhi,  P.  Givi  and  S.B.  Pope,  “Implementation  of  the  Velocity-Scalar 
Filtered  Mass  Density  Function  for  Large  Eddy  Simulation  of  Turbulent 
Reacting  Flows,”  Bulletin  of  the  American  Physical  Society ,  50,  p.  261,  58th 
Annual  Meeting  of  the  Division  of  Fluid  Dynamics  of  the  American  Physical 
Society,  Chicago,  IL,  November  20-22,  2005. 

2.  T.G.  Drozda,  M.R.H.  Sheikhi,  P.  Givi  and  S.B.  Pope,  “Large  Eddy  Simulations 
of  a  Bluff-Body  Stabilized  Hydrogen-Methane  Jet  Flame’  Bulletin  of  the 
American  Physical  Society ,  50,  p.  101,  58th  Annual  Meeting  of  the  Division 
of  Fluid  Dynamics  of  the  American  Physical  Society,  Chicago,  IL,  November 
20-22,  2005. 

3.  Invited:  P.  Givi,  “Quality  Assessment  of  the  Filtered  Density  Function  for 
Large  Eddy  Simulation,”  Proceedings  of  the  First  Workshop  on  Quality 
Assessment  of  Unsteady  Methods  for  Turbulent  Combustion,  Available  on 
CD  only,  Darmstadt,  Germany,  June  16-17,  2005. 

4.  M.R.H.  Sheikhi,  P.  Givi,  and  S.B.  Pope,  “Latest  Developments  in  FDF 
Formulation,”  Proceedings  of  the  Fourth  Joint  Meeting  of  the  U.S.  Sections 
of  the  Combustion  Institute,  available  on  CD  only.  Philadelphia,  PA,  March 
20-23,  2005. 
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5.  M.R.H.  Sheikhi,  P.  Givi  and  S.B.  Pope,  “Joint  Velocity-Scalar  Filtered  Mass 
Density  Function  for  Large  Eddy  Simulation  of  Turbulent  Reacting  Flows,” 
Bulletin  of  the  American  Physical  Society,  49(9),  p.  217,  57th  Annual  Meeting 
of  the  Division  of  Fluid  Dynamics  of  the  American  Physical  Society,  Seattle, 
WA,  November  21-23,  2004. 

6.  P.  Givi,  M.R.H.  Sheikhi,  T.G.  Drozda,  and  S.L.  Yilmaz,  “A  New  Strategy 
for  Turbulence  Simulation,”  Special  Session  on  Multiscale  Algorithms 
in  Computational  Fluid  Dynamics,  Meeting  #1002  of  the  American 
Mathematical  Society,  Pittsburgh,  PA,  November  6-7,  2004. 

7.  T.G.  Drozda,  M.R.H.  Sheikhi,  P.  Givi,  F.A.  Jaberi  and  S.B.  Pope,  “Scalar 
Filtered  Density  Function  for  Large  Eddy  Simulation  of  a  Piloted  Diffusion 
Flame,”  Session  HT-10E,  2003  ASME  International  Mechanical  Engineering 
Congress,  Washington,  D.C.,  November  15-21,  2003. 

8.  T.G.  Drozda,  M.R.H.  Sheikhi,  P.  Givi,  and  S.B.  Pope,  “Large  Eddy  Simulation 
of  Scalar  Transport  in  a  Turbulent  Mixing  Layer,”  Proceedings  of  the  2003 
Fall  Technical  Meeting  of  the  Combustion  Institute,  Eastern  States  Section, 
in  Chemical  and  Physical  Processes  in  Combustion,  pp.  237-240,  University 
Park,  PA,  October  26-29,  2003. 

9.  M.R.H.  Sheikhi,  T.G.  Drozda,  P.  Givi,  F.A.  Jaberi,  and  S.B.  Pope,  “Large 
Eddy  Simulation  of  Sandia  Flame  D,”  Proceedings  of  the  2003  Fall  Technical 
Meeting  of  the  Combustion  Institute,  Eastern  States  Section,  in  Chemical  and 
Physical  Processes  in  Combustion,  pp.  241-244,  University  Park,  PA,  October 
26-29,  2003. 


Invited  Public  Lectures  and  Seminars: 


During  the  course  of  this  research,  the  PI  delivered  the  following  invited  lectures 
and  seminars. 

1.  Affordable  and  Reliable  Prediction  of  Turbulent  Combustion,  Department  of 
Mechanical  Engineering,  University  of  Utah,  Salt  Lake  City,  UT,  February  3, 
2006. 

2.  Affordable  and  Reliable  Prediction  of  Turbulent  Combustion,  Computational 
Aerosciences  Branch,  NASA  Langley  Research  Center,  Hampton,  VA, 
December  15,  2005. 


21 


3.  Affordable  and  Reliable  Prediction  of  Turbulent  Combustion,  The  Center  for 
21st  Century  Energy,  FALL  2005  Seminar  Series,  Department  of  Mechanical 
Engineering,  Massachusetts  Institute  of  Technology,  Cambridge,  MA,  October 
14,  2005. 

4.  Introduction  to  Simulation  and  Computational  Fluid  Dynamics,  16th 
International  Activated  Carbon  Conference  and  Courses,  Pittsburgh,  PA, 
October  6,  2005. 

5.  LES  of  High  Speed  Combustion,  The  Media  Effects  Project,  Presentation  for 
Office  of  Secretary  of  Defense,  NASA  Langley  Research  Center,  Hampton,  VA, 
October  4,  2005. 

6.  Subgrid  Scale  Modeling  for  Large  Eddy  Simulation  of  Turbulent  Reacting 
Flows:  A  Review,  Institute  for  Scientific  Computing,  University  of  Wyoming, 
Laramie,  WY,  September  15,  2005. 

7.  Subgrid  Scale  Closure  in  Turbulent  Combustion:  A  Review,  Wright  Patterson 
Air  Force  Base,  AFRL/VA,  Dayton,  OH,  August  24,  2005. 

8.  Filtered  Mass  Density  Function  for  Subgrid  Scale  Modeling  of  Turbulent 
Diffusion  Flames,  ARO/AFOSR  Contractors’  Meeting  in  Chemical 
Propulsion,  Indianapolis,  IN,  June  22,  2005. 

9.  Quality  Assessment  of  the  Filtered  Density  Function  for  Large  Eddy 
Simulation,  First  Workshop  on  Quality  Assessment  of  Unsteady  Methods  for 
Turbulent  Combustion,  Darmstadt,  Germany,  June  17,  2005. 

10.  Subgrid  Scale  Closure  in  Turbulent  Combustion:  A  Review,  CNRS,  Ecole 
Centrale,  Paris,  France,  June  13,  2005. 

11.  Subgrid  Scale  Closure  in  Turbulent  Combustion:  A  Review,  Colloquium  of 
Sibley  School  of  Mechanical  and  Aerospace  Engineering,”  Cornell  University, 
Ithaca,  NY,  March  29,  2005. 

12.  FDF  and  SPARK  for  LES,  The  Media  Effects  Project,  Presentation  for  Office 
of  Secretary  of  Defense,  NASA  Langley  Research  Center,  Hampton,  VA,  March 
11,  2005. 

13.  Mathematical  Modeling  and  Computer  Simulation  of  Combustion  and 
Propulsion,  Poster  presented  at  the  Energy  Event:  A  Pitt  Showcase  of  Energy 
Innovation,  University  of  Pittsburgh,  PA,  December  9,  2004. 
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14.  Implementation  of  FDF  in  SPARK-3D,  The  Media  Effects  Project, 
Presentation  for  Office  of  Secretary  of  Defense,  NASA  Langley  Research 
Center,  Hampton,  VA,  November  2,  2004. 

15.  Subgrid  Scale  Closure  for  Large  Eddy  Simulation  of  Turbulent  Combustion, 
Advanced  Analysis  Processes,  Rocketdyne  Propulsion  &  Power,  The  Boeing 
Company,  Canoga  Park,  CA,  July  22,  2004. 

16.  Fire  and  Its  Mystique:  From  the  Caveman  to  the  Rocket  Scientist,  Inaugural 
Lecture  for  Endowed  WILLIAM  KEPLER  WHITEFORD  Professorship, 
University  of  Pittsburgh,  PA,  February  19,  2004. 

17.  A  New  CFD  Strategy  for  Simulation  of  Turbulence  and  Combustion,  Center 
for  Bioelectronics,  Biosensors  and  Biochips,  Department  of  Mechanical 
Engineering,  Virginia  Commonwealth  University,  Richmond,  VA,  December 
10,  2003. 

18.  Supercomputer  Simulations  in  Combustion  and  Propulsion:  Making  it 
Much  Easier  to  be  a  Rocket  Scientist,  Poster  presented  at  the  Technology 
Commercialization  Alliance,  University  of  Pittsburgh,  PA,  December  9,  2003. 

19.  Subgrid  Scale  Modeling  in  Turbulent  Combustion:  A  Review,  39th 
AIAA/ASME/SAE/ASEE  Joint  Propulsion  Conference^  Exhibit,  Huntsville, 
AL,  July  23,  2003. 

20.  Filtered  Mass  Density  Function  for  Subgrid  Scale  Modeling  of  Turbulent 
Diffusion  Flames,  ARO/AFOSR  Contractors’  Meeting  in  Chemical 
Propulsion,  Williamsburg,  VA,  June  25,  2003. 

21.  A  Novel  CFD  Strategy  for  Turbulence  and  Combustion  Simulation, 
Department  of  Mechanical  and  Industrial  Engineering,  University  of  Toronto, 
Toronto,  Ontario,  Canada,  March  21,  2003. 

22.  Large  Eddy  Simulation  of  Turbulent  Combustion,  National  Energy  Technology 
Laboratory,  U.S.  Department  of  Energy,  Morgantown,  WV,  February  11,  2003. 

23.  A  Novel  CFD  Strategy  for  Turbulence  Simulation,  Institute  for  Complex 
Engineering  Systems  and  Department  of  Mechanical  Engineering,  CMU, 
Pittsburgh,  PA,  January  24,  2003. 
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(c) 


Figure  1:  Contour  surface  of  the  (4>) L  field  in  the  3D  mixing  layer  at  t  =  80  as 
obtained  by:  (a)  DNS,  (b)  VSFMDF  and  (e)  Srnagorinsky. 
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Figure  2:  Cross-stream  variations  of  some  of  the  Reynolds-averaged  components  of 
r/,  at  t  =  80.  The  thick  solid  and  thin  dashed  lines  denote  LES  predictions  using 
VSFMDF  and  Smagorinsky  closures,  respectively.  The  circles  show  the  filtered  DNS 
data. 
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Figure  3:  Cross-stream  variations  of  some  of  the  components  of  R  at  t  =  80. 
The  thick  solid  and  thin  dashed  lines  denote  LES  predictions  using  VSFMDF  and 
Smagorinsky  closures,  respectively.  The  circles  show  the  filtered  DNS  data. 
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Figure  4:  Cross-stream  variations  of  r  at  t  =  80.  The  thick  solid  and  thin  dashed 
lines  denote  LES  predictions  using  VSFMDF  and  Smagorinsky  closures,  respectively. 
The  white  and  black  circles  show  the  filtered  and  unfiltered  DNS  data,  respectively. 
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Figure  6:  Time  averaged  recirculation  features  as  predicted  by  experiments  and 
LES.  The  color  contours  denote  the  streamwise  velocity.  Superimposed  are  the 
streamlines  and  the  velocity  vectors. 
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Figure  7:  Time  averaged  recirculation  features  as  predicted  by  LES.  The  color 
contours  describe  the  strearnwise  velocity.  Superimposed  are  the  streamlines  and 
the  velocity  vectors. 
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(a)  LES-50 


(b)  LES-8-5 


(c)  LES-143 
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Figure  8:  Instantaneous  (a,  b,  c),  and  time  averaged  (d,  e,  f),  contours  of  the  mixture 
fraction  as  predicted  by  LES  of  several  non-reacting  flow  configurations.  The  six 
iso-surfaces  arc  0.06  (blue),  0.13,  0.34,  0.45,  0.63,  0.90  (white).  The  plot’s  axes  are 
normalized  by  the  diameter  of  the  fuel  stream. 
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(a)  HMl 


(c)  HM1E-S2 


Figure  9:  Time  averaged  recirculation  features  as  predicted  by  SFMDF.  The 
contours  denote  the  streamwise  velocity.  Superimposed  are  the  streamlines  and 
the  velocity  vectors. 
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Figure  10:  Radial  profiles  of  the  mean  value  of  the  streamwise  velocity  for  the 
reacting  How.  - SFMDF,  •  HM1,  ■  HM1E-S1,  ♦  HM1E-S2. 
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7<  z  >  7<  z  >  7<  z  > 


Figure  11:  Radial  profiles  of  the  mean  values  of  the  mixture  fraction  for  the  reacting 
flow.  - SFMDF,  •  HM1. 
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Figure  12:  Radial  profiles  of  the  mean  values  of  the  temperature  for  the  reacting 
flow.  - SFMDF,  •  HM1. 
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Figure  13:  Radial  profiles  of  the  RMS  values  of  the  temperature  for  the  reacting 
flow.  - Resolved  SFMDF,  -  -  Total  SFMDF,  •  HM1. 
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Figure  14:  PDFs  of  the  resolved  mixture  fraction  at  (top- to- bottom  rows)  x/Db  = 

0.3,  0.6,  0.9,  1.3,  1.8  and  selected  radial  locations  (r  [mm]).  -  SFMDF, 

HM1. 
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Figure  15:  PDFs  of  the  resolved  temperature  at  (top-to-bottom  rows)  x/Db  = 

0.3,  0.6,  0.9,  1.3,  1.8  and  selected  radial  locations  (r  [mm]).  -  SFMDF, 

HM1. 
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